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ABSTRACT 

We show that the hierarchical model at finite volume has a symmetry group 
which can be decomposed into rotations and translations as the familiar Poincare 
groups. Using these symmetries, we show that the intricate sums appearing in 
the calculation of the high-temperature expansion of the magnetic susceptibility 
can be performed, at least up to the fourth order, using elementary algebraic 
manipulations which can be implemented with a computer. These symmetries 
appear more clearly if we use the 2-adic fractions to label the sites. We then apply 
the new algebraic methods to the calculation of quantities having a random walk 
interpretation. In particular, we show that the probability of returning at the 
starting point after m steps has poles at D = —2, —4, .... — 2m , where D is a free 
parameter playing a role similar to the dimensionality in nearest neighbor models. 



1. Introduction 



The hierarchical model played an important role in the development of the 
ideas of the renormalization group f 1 because the renormalization group transforma- 
tion for this model reduces to a simple integral equation which has been studied 
in great detail.' 31 In a recent paper) 41 we calculated numerically the magnetic sus- 
ceptibility of the hierarchical Ising model as a function of the temperature and for 
two values the parameter used in the e-expansion. Calculations have been carried 
with up to 2 18 sites. We found that the numerical data can be fitted precisely with 
a simple power law of the form (1 — f3/ 13q)~ 9 in the whole high-temperature region, 
i.e for (5 G [0, f3 c ]. As a consequence, it is possible to obtain a simple approxi- 
mate formula for the high-temperature coefficients of the susceptibility in terms 
of g and (3q. This approximate formula implies approximate relations 14 ' 51 among 
the coefficients which deserve a systematic investigation. Unfortunately, due to 
the non-locality of the hierarchical model, the calculation of the high-temperature 
expansion is quite non-trivial even for the very first orders. 

The goal of this article is to explain in a pedestrian way how the symmetries of 
the hierarchical model can be used to reduce the convoluted sums appearing in the 
high temperature expansion, at least up to order 4 in /3, to simple algebraic manip- 
ulations which can be easily implemented with a computer. In section 2, we show 
that even at finite volume, the hamiltonian of the hierarchical model is invariant 
under a group of transformation which can be decomposed into translations and 
rotations as the familiar Poincare groups. 

One important reason to work at finite volume is that high-temperature co- 
efficients of arbitrary order can be calculated with an arbitrary precision using a 
computer, if the number of sites is not too large. This can be done by introducing 
a f3 expansion within the exact method of integration described in Ref. [4]. Since 
it is very easy to make combinatoric mistakes in calculating high-temperature co- 
efficients, formulas valid for arbitrary volume allow unambiguous checks if we can 
compare them with sufficiently precise numerical results obtained with a finite 
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number of site. A numerical program performing this calculation is provided in 
Appendix 2. This program can be used to check the validity of the analytic results 
presented hereafter with at least 10 significant figures, when the number of sites 
does not exceed 2 9 . More fundamentally the renormalization group approach re- 
quires techniques of integration at finite volume and the results obtained here can 
be used in this context. 

In section 3 , we use extensively the translation invariance to calculate the 
first four coefficients of the magnetic susceptibility. The final result is that each 
of these coefficients can be written a sequence of products and Fourier transforms. 
In section 4, we show that these sequences of operations can be reduced to simple 
algebraic manipulations which can be implemented on a computer. A ready-to-use 
Mathematica program is provided in Appendix 1. The methods presented here can 
be used for higher order calculations, however it is not clear that the calculation 
of coefficients of an arbitrary order can be completely reduced to simple algebraic 
manipulations. This issue is briefly discussed in the conclusions. 

The results of section 2 and 4 could have been derived more elegantly if we had 
labeled the sites using 2-adic fractions as in Ref. [7]. This is explained in section 
5 which, unlike section 2, assumes some familiarity with the mathematics involved 
in this reformulation. However, section 5 is not essential to the understanding of 
the other sections. 

The results obtained in section 4 take a more compact form in the infinite 
volume limit. In section 6, we give the explicit form of each terms appearing 
in the calculation of the first four coefficients in this limit. These expressions 
have poles at unphysical values of D, a parameter controlling the strength of the 
bilinear couplings among the spins (see section 2), and playing a role similar to the 
dimension in nearest neighbor models. By unphysical values of D, we mean values 
such that changing the sign of one spin in one of the configuration where all the 
spins are aligned causes an infinite increase in energy. Incidently, all the quantities 
calculated in that section have a random walk interpretation, provided that a 
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proper normalization is used. In general, the calculation the m-th coefficient of the 
high temperature expansion involves the calculation of the probability of returning 
at the starting point after m steps and the probability of visiting less than two 
sites in m steps. In the infinite volume limit, these quantities can be calculated 
for arbitrary m. Their poles are located at D — —2, —4, —2m and D = — 
respectively. 

2. The Hierarchical Ising Model and its Poincare Group 

In this section we describe the hierarchical Ising model and its invariance under 
a group of transformation which can be decomposed into translations and rotations 
in a way similar to the Poincare group used in relativistic quantum field theory. 
Hierarchical models' 11 are specified by a non-local hamiltonian bilinear in the spin 
variables and a local measure of integration which will not be specified in this 
section. The specific choice of an Ising measure where the spins take only the 
values ±1 will only be made in the next sections. We first recall the form of the 
non-local hamiltonian. 

The hierarchical models considered here require the number of sites to be 2 n . 

For convenience, we label the sites with n indices x n xi, each index being or 1 

. In order to give a concrete meaning to this notation, one can divide the 2 n sites 
into two boxes, each containing 2 n_1 sites. If x n = 0, the site is the first box, if 
x n — 1, the site is in the second box. Repeating this procedure n times (for the two 
boxes, their respective two sub-boxes, etc.), we obtain an unambiguous labeling for 
each of the sites. In the following, we often use the symbol i as a short notation 
for the sequence x n x\. 

The hamiltonian of the hierarchical model reads 
1 n 

h =-2^ 1 e ( e ^ M )) 2 - (2- 1 ) 

1=1 X n ,...,Xl+l Xl,....,Xl 
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The model has a free parameter c for which we shall use the parametrization 



c = 2 1 ~d. 



(2.2) 



The parameter D plays a role similar to the dimension in nearest neighbor models. 
The parameter of the epsilon-expansion can be defined as e = 4 — D. When 
D > 4, the model has a trivial continuum limit. When D < 2, the model does 
not have a phase transition at finite temperature.' 11 These two rigorous results can 
be understood heuristically in terms of the the self-intersection properties of the 
random walk associated with H, by noticing that the Haussdorff dimension of this 
random walk is 2/D. 

For a given / in Eq. (2.1), the first sum can be interpreted as a sum over boxes 
of size 2 l , while the second sum is over the spins inside each of these boxes. Given 
two distinct sites x and y, the corresponding spins have an interaction whenever 
I is large enough to have x and y within the same box. It is clear that if x 
and y belong to the same box of size 2 l , they also belong to the same boxes of 
larger size, until the maximal size is reached. Since the respective contributions 
to the interaction are (|)' in Eq. (2.1), the total interaction can be expressed as 
a truncated geometrical series. In order to describe quantitatively this situation, 
we define a function v(x,y) which indicates the "level" / at which x and y start 
to differ. More precisely, if x and y are distinct, v(x,y) = I when x m = y m for 
all m such that n > m > I > and x\ ^ y\. At coinciding arguments, we define 
v(x,x) = 0. The hamiltonian function H can then be rewritten as 




(2.3) 



where 




(2.4) 
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and 



L = ((|)-(^r +1 )(i-|) 



(2.5) 



As made clear by the above equations, the strength of the interaction between 
two spins a x and a y depends only on the value of v(x,y). Consequently the 
invariance of v(x,y) under a group of transformation (see Eq. (2.10)) implies 
the invariance of H under corresponding transformations, in a way which will be 
made precise at the end of this section. The invariance of v(x, y) under a group of 
transformation appears rather clearly in the reformulation of the hierarchical model 
proposed in Ref . [7] , where the sites are associated with 2-adic fractions which can 
be transformed into each other using addition and multiplication. However, it is 
possible to obtain the results that will be needed later using only a few basic results 
in arithmetic modulo 2 n . The reader familiar with the 2-adic numbers will notice 
that since the 2-adic integers can be defined 1101 as the projective limit of the integers 
modulo 2 n , the finite volume presentation given here is implicit in the formulation 



In order to describe the invariance of v(x,y), we shall associate with the se- 
quence of 0's and l's x n x\, a rational number of the form 



Since this defines a one-to-one correspondence, we shall use the same symbol x for 
the sequence of and 1 and the rational number. We can now define additive and 
multiplicative operations which will be used to obtain a group of transformation. 

If two numbers x and y have the form given in Eq. (2.6), then x + y can also 
be written as a rational number where the denominator is 2 m with m < n. If we 
drop the integer part of x + y, we obtain a rational number having the form of Eq. 
(2.6). Equivalently, we can write x = q/2 n and y = r/2 n with q and r integers 
between and 2 n — 1 and add q and r modulo 2 n . Since the integer modulo 2 n 



of Ref. [7]. 



n 




(2.6) 



m=l 
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form an additive group, the set of fractions associated with the sites form a group 
for the addition modulo 1. In order to clarify the ideas, let us give an example: 

the additive inverse of 2~ n modulo 1 can be written as 2~~ n + 2~ n+1 + + 2 _1 . 

This makes clear that the addition modulo 1 is distinct from the addition obtained 
by adding each of the indices (x m + y m ) modulo 2. Furthermore, the odd integers 
modulo 2 n form a multiplicative group. We can pick a canonical form for the 
representatives of such integers as 

u = 1 + 2z (2.7) 

where z is a positive integer between and 2 n_1 — 1. Obviously, if x has the form 
(2.6), then ux has also the form (2.6) after discarding its integer part. 

We are now in position to define a group of transformation acting on the 
fractions associated with the sites. If x and a have the form of Eq. (2.6) and u has 
the form of Eq. (2.7), we define a transformation of x depending on a and u and 
denoted x[u, a] which reads 

x[u, a] = ux + a , (2.8) 

where the r.h.s is understood modulo 1. It is clear that u — I, a — gives the 
identity transformation. From the composition law (x[u, a])[v, b] = x[uv,va + b), 
one sees that an inverse transformation is obtained for v = u" 1 and b = —u^ 1 a. 
Consequently these transformations form a (non-abelian) group. 

We can interpret x[0, a] as a translation and x[u, 0] as a rotation like in the usual 
Poincare groups. In that sense, this is a "global" group of transformation. This is 
in contrast with the symmetries noted by Dyson' 11 which consists in interchanging 
x m xi + ilxi-i xi and x m z/ + iCte/_i x\ "locally". 

In the following, we shall frequently use the fact that x[u, a] with u and a 
fixed is a one-to-one map. The proof of this statement is simple. If ux + a = 
uy + a modulo 1, then u(x — y) = modulo 1, writing x — y as q/2 n , we obtain 
that q = modulo 2 n and consequently, x = y modulo 1. 
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We can now make a connection between v(x, y) and x — y through the following 
result: 



v(x,y)—l x — y = q/2 modulo 1 for some odd integer q (2-9). 

In order to prove this, we write x = r/2 n and y = s/2 n with r and s integers 
between and 2 n — 1 and we see that r.h.s of the equivalence relation means that r 
and s are equal modulo 2 n ~ l but differ modulo 2 n ~ l+1 . Returning to the canonical 
form of Eq. (2.6), we see that is equivalent to v(x,y) = I. Noticing that x — y is 
translation invariant and that the product of odd integers is still an odd integer, 
we obtain the invariance of v(x, y) under the transformations (2.8) in the following 
way: 

v(x[u,a],y[u,a]) =v(x,y) . (2.10) 
Since K. xy depends only on v(x,y), this immediately implies that 

K^y = ^- X [u,a]y[u,a] ■ (2-H) 

It is time to state the main result of this section: The hamiltonian H defined 
by Eq. (2.1) is invariant under the transformation 

°x > (*x[u,a] ( 2 - 12 ) 

for any a as in Eq. (2.6) and any u as in Eq. (2.7). 

This follows from the fact that x[u, a] defines one-to-one map, consequently we can 
rewrite the transformed hamiltonian in the following way 

y ] ^■xyO'x\u,a\ (T y\uM = ^ 1 ^■x[u- 1 ,-u- 1 a]y[u- 1 ,-u- 1 a\ <T xO'y ( 2 -13) 

x,y x,y 

and then use the invariance of K xy shown in Eq. (2.11). 
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As is well-known, the invariance of H under a group of transformations implies 
identities for the correlation functions. In order to define these, we need to intro- 
duce a local measure of integration. In the following sections, we shall consider the 
case of an Ising measure and define 

Z= e ^ H ( 2 - 14 ) 

{a=±l} 

and 

< <j x a y >= Z- 1 Wy^ 1311 , (2-15) 

{a=±l} 

We can then prove that 

< °x[u,a]°y[u,a] > = < O x O y > (2.16) 

This follows from the fact that we can always relabel the integration variables 
according to Eq. (2.12.) and then use the invariance of H. This argument is 
clearly independent of the choice of the local measure. 

3. The High- Temperature Expansion 
of the Hierarchical Ising Model 

We now proceed to the calculation of the magnetic susceptibility per site for 
the hierarchical Ising model. We define the magnetic susceptibility per site as 



Xr 



X 

Using translation invariance, we find that 

XnW) = < <7o(X>*) > > ( 3 - 2 ) 

x 

the choice of being arbitrary. In the case of the Ising model we can use the fact 
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that <7q = 1 and write 

X n(P) = 1 + < a0CTx > ' ( 3 - 3 ) 

x;x^O 

Using the high-temperature expansion of < (Tq(t x >, we can calculate the 
coefficients of the expansion 

Xn (P) = 1 + h,nP + b 2 , n P 2 + - (3.4) 

A brief review of the standard methods used to perform such a calculation with 
arbitrary long-range Ising models is given at the beginning of section II of Ref. [6] . 
For the sake of briefness, we shall only recall the final result: < <jq<j x > can be 
written as the ratio of two expansions in tanh((3K. yz ). The terms of each of them 
correspond to graphs where a given link y — z can only appear once (in which case 
we have a have a factor tanh(j3~K yz )) or zero times (in which case we have a factor 
1). At the denominator, all the sites must be visited an even number of times. At 
the numerator, all the sites must be visited an even number of times except for 
and x which are visited an odd number of times. This provide an expansion in 
tanh(f3~K yz ) which will require further expansion in order to be recast in the form 
(3.4). It has to be noted that for hierarchical models, the use of Feynman diagrams 
in the recursion formula can be reduced to a purely combinatoric problem (no 
sums or integrals), while the high-temperature expansion requires the evaluation 
of nested sums - a situation which somehow the opposite of what we encounter for 
nearest neighbor models. We now proceed order by order to the calculation of the 
b m ,n for m = 1, 2, 3, 4 and n arbitrary. 

The first order contributions to < uqu x > comes from high-temperature 
graphs with one link joining to x. The contribution of these graphs is 

J2^nh(PK 0x ) . (3.5) 

X 

Note that K. xy is zero when x = y and we do not need to subtract the contribution 
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with x — 0. We thus obtain 



h,n = V Koa; . 



(3.6) 



X 



Proceeding similarly with graphs with two links, we obtain 



yx 



(3.7) 



y,x 



At order three, we need to consider the (3 3 terms of Eq. (3.5) in addition to the 
lowest order contributions from the graphs with three links. The result is 



Up to now, we did not have to take into account the contributions from the denom- 
inator because the first non-trivial contributions to Z come from triangles which at 
lowest order give order (3 3 terms. Since < ao<J x > yields contributions which are 
at least of order (3 when x ^ 0, the denominator has to be taken into account at 
order 4. These contributions cancel the lowest order of the graphs with a triangle 
at one or the other end and all which is left are minus the graphs where the link 
— x is also occupied by one side of a triangle (since such a graphs do not appear 
at the numerator). Adding to this the lowest order contribution from the graphs 
joining 5 distinct sites and the order /5 4 of the graphs with two links, we obtain 






y ,z,w,x 
all distinct 



xy 



(3.9) 



3 5^( K 0y) 3 K yx . — -^2K y(K yx ) 3 + 77 y^(K 0:r ) 4 




Note that we have taken into account that K T7 , = K. 



yx- 
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From the translation invariance of K. xy , it is clear that K. xy is a function of 
x — y only. Consequently, we can write 

2 n -l 

K *v = Yn £ KC^e-^-y) . (3.10) 

fc=0 

The explicit form of K(fc) will be given in Eq. (4.12). Our task in now to rewrite 
the expressions for the b m ,n obtained in the previous section as unrestricted sums 
over y, z, ...x with appropriate compensating terms, plug the expansion (4.1) in 
them and use the well-known results 

2 n -l 

J2* XZ G zy = - £ K(k)G(k)e-^ k ^ . (3.11) 

z k=0 

and 

2 n — 1 

_L ^ e -m q = Sko (312) 

Proceeding this way, we obtain 

6i,n = K(0) (3.13) 



6 2 , n = (K(0)) 2 -K 2 (0) (3.14) 
6 3 , n = (K(0)) 3 - 2K(0)K2(0) + ^K3(0) - (R3)(0) (3.15) 

6 4 , n =(K(0)) 4 + ^K(0)K3( ) - 1k4( ) + 2(K2K2)(0) 

- 3R2(0)(K(0)) 2 + 3(K2(0)) 2 - £ 4 (0) - 2K(0)((lo 3 )(0) 

Note that expressions of the same form but with different coefficients would be 
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obtained if instead of using a Ising measure, we had an arbitrary measure. This 
can be seen in general by using the random path representation.' 111 



4. Explicit Expressions for the High- Temperature Coefficients 

Our next task consists in calculating explicitly the above expressions. For this 
purpose, we first introduce the function N(k) defined over the integers modulo 2 n : 

„ r ,, s f 2~ l if k can be divided by 2 l but not by 2 l+1 , t 

N(k) = < j j 

l2- n if k = 0modulo2 n 

This function will appear in K.(k). We now introduce short notations for some 
functions of x 

2 -av{x,0) — > x g 

S x ,o — > S (4.2) 



and for some functions of k 



In addition, we also define 



{N{k))o s 

ho — A (4.3) 



P(l + a) = ^T-^ (4.4) 

v 1 1 - 2~( 1+a ) 



G ( i + a )= 1 _ 1 2 X ffl ) (4 - 5) 



Q(l+a) = 2-( n+1 )( a+1 ) (4.6) 
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With these notations, 

K x0 — > P{\ + 1)(X 1+ * - Q(l + ^)* - (1 - Q(l + ^))<S) (4.7) 

It has to be noted that Q(l + a) vanishes in the infinite volume limit. Taking this 
limit simplifies the calculation, but obviously prevents us from checking it using 
exact results at finite volume. 

We can now manipulate the symbols introduced above by just using their mul- 
tiplication table and the rules to perform their Fourier transform. The multiplica- 
tions appearing in Eqs. (3.13-16) can be performed as a multiplication denoted * 
(respectively * in Fourier transform) in an abstract associative algebra where the 
basis of the underlying vector space is given by the symbols introduced in Eq. (4.2) 
(resp. Eq. (4.3)). Since in these two equations s can be understood as a continuous 
index, this vector space is infinite dimensional. However, for calculations at finite 
order, only finite dimensional subspaces are used. The multiplication table for the 
X s and 5 is 

x s * X t = X s+ t 

X s *5 = 5 (4.8) 
5*5 = 5 

Similarly for K. s and A, we have 

fC s *fCt = fC s+ t 

K a * A = 2' ns ^A (4.9) 
A * A = A 

A straightforward but tedious calculation yields the Fourier transforms 

Xs = (Gis))- 1 ^^ - 2 S - 1 /C ) + 2 n Q(s)A 

Xq = 2 n A (4.10) 

5 = 1 
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and their inverse 

K a = G{\ + a){X 1+a - Q(l + a)Xo) + 2 a 5 

K Q = 5 (4.11) 

A = 2~ n X Q 

The basic ingredient in this calculation is the identity (3.12). However, a more 
organized calculation can be performed using general results given in Ref. [9] (see 
next section). With these results, one can calculate K(/c) appearing in Eq. (3.10) 

K(fc) =P(l + j 5 )(G(l + ±))- 1 . 

(JCi + (~P(1 + ^) + Q(l + ±)G(1 + ^))/Co) 

This results can also be used for Feynman diagram calculations. 

Finally, we need to evaluate all these functions when their argument is 0. This 
again can be accomplished in a purely algebraic fashion, by multiplying by A (resp. 
S) when an odd (resp. even ) number of Fourier transform has been performed and 
retaining the coefficient of A (resp. 5). 

Using these results we find (recalling the definition of c given in (2.2)) 

h, n = a - lr\\{i - (|d(i - ^r 1 - (2 n - i){\r +l ) (4.i2) 

and 

b2, n =(w 2 - a - f r 2 [(f ) 2 a - 4 )n)(1 - ¥ rl 

- 2(^+ 2 (l - (£)»)(! - £)"! - (2- - 1)(^ +1 )] 

The next two coefficients are more involved and writing them explicitly might 
not be the most efficient way to communicate the information to the reader. In- 
stead, we give in Appendix 1 a Mathematica program which calculates the first 
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(4.12) 



(4.13) 



four coefficients and allows further symbolic or numerical manipulations. This pro- 
gram has been tested in various ways. The most convincing check is the precise 
comparison between the numerical values obtained by replacing D and n in the 
expression of the coefficients and the exact values obtained by introducing a high 
temperature expansion in the recursive integration used in Ref. [4]. A program 
performing this calculation is given in Appendix 2. 

Before closing this section, let us recall that in order to obtain a physically 
interesting model) 11 it is essential to keep bi :U finite in the limit where n becomes 
infinite. Otherwise, flipping one spin in one of the configuration where all the spins 
are aligned causes an infinite change in energy, this requirement imposes \c\ < 2. If 
in addition we require ferromagnetic couplings, we see from Eq. (2.4) that c must 
be real, strictly positive and less than four. Combining the two requirements yields 
< c < 2 which incidently corresponds to D real, strictly positive and finite. 



5. Remarks Concerning the 2-adic 
Formulation of the Hierarchical Model 

Some of the results presented in the previous sections can derived more ele- 
gantly if we had started immediately with the reformulation of Ref. [7]. In this 
section, we briefly explain how this can be done. Unlike the previous sections, this 
one assumes some familiarity with the p-adic numbers. 

The function 2 V ^ X,V ^> introduced in section 2, is a regularized version of the 
2-adic distance. Namely, 



, \x — yh if \x — yh > 1 
2 v(x,y) — ^ yiz 1 yiz (5.!) 

1 if \x — 2/j 2 < 1 



Using this, the invariance of v(x,y) given in Eq. (2.10) follows easily. 

Similarly, N(k) introduced in section 4, is a regularized version of the 2-adic 



1G 



norm restricted to the 2-adic integers. 



N(k) 



\k\ 2 if 1 > \k\ 2 > 2 
2~ n if |jfe| 2 < 2- n 



—n 



(5.2) 



Given this, one can use the results of Ref. [9] to calculate the Fourier transforms. 
In particular, the function G(l + a) is the gamma function associated to a set of 
multiplicative characters. 

6. The Poles in the c-Complex Plane, in the Infinite Volume Limit 

In the infinite volume limit, the expressions for the coefficients simplify signif- 
icantly. This limit will always be understood in this section and we shall omit the 
reference to n. In taking this limit, we assume < c < 2. For instance, for the 
first coefficient, we obtain 



From Eq. (2.4), it is obvious that for any tti, bm will have a zero of order ui at 
c = and a pole of order m at c = 4. Also, from the procedure used to calculate 
the coefficients, one can see that b m will have a (&i) m term (which would be the 
only contribution in the gaussian case). This term has a pole of order m at c = 2. 
In order to get rid of these zeroes and poles, we shall consider b m /(b™) in the 
calculations below. In addition, with this normalization, the quantities calculated 
have a random walk interpretation when < c < 2. This is due to the fact that 
K-cy/ J2 Z Kq z can then be seen as the probability for going from x to y in one step. 

A quantity which is easy to calculate is the probability for visiting no more 
than two sites, including the starting point, after m step. This quantity reads 
2~2 x (K-0x) m /(2~2z Koz) m and appears in the calculation of b m /(bi) m . After reex- 
pressing it in terms of the Fourier transform and using the techniques developed 




X 



(6.1) 
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in section 4, we obtain 

K^(0) (2-c) m 



(K(0)) m 2- 1+2m -c m 

Proceeding similarly for the other quantities entering in &3/(&i) 3 and 64/ (61) 4 , 
we obtain 

K 3 (0) 6(-2 + c) 3 c 2 



(K(0)) 3 (8-c 2 ) (-16 + c 3 )' 
K 4 (0) (-2 + c) 4 (-256 - 160 c 2 - 14 c 5 + 16 c 3 ) 



and 



(K(0)) 4 ("8 + c 2 ) (-16 + c 3 ) (-32 + c 4 ) 



;iT 2 K 2 ) (0) (-2 + c) 4 (-128 c 2 - 64 c 3 + 12 c 5 ) 



(6.3) 



(6.4) 



kT 4 (0) (-8 + c 2 ) 2 (8 + c 2 ) (-32 + c 3 ) 



(6.5) 



Note that Eqs. (6.3) and (6.4) can be generalized to the calculation of the 
probability for a return at the starting point after m steps. A detailed calculation 
shows that this quantity can be expressed as 



K™(0) 
(K(0)) m 



l—n \ / 



This shows that unless unexpected cancelation occurs, b m will have poles at c = 

2 1+1 , 2 1+1 / 2 , 2 1+1 /™ or in other words at D = -2, -4, , -2m. On the other 

hand, the quantity appearing in Eq. (6.6) has a zero of order m at c = 2. This is 
not obvious from the r.h.s, however it can be seen by expanding the denominator 
about 2 and summing over /, that this is the case for each term of the resulting 
expression. The location of these poles plays an important role if one tries M to use 
the 1/D-expansion to calculate the corrections to the gaussian contribution. 
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7. Conclusions 



We have shown that the symmetries of the hierarchical model can used exten- 
sively to calculate the high temperature expansion of the magnetic susceptibility. 
Up to order 4, the calculation can be performed using algebraic methods which can 
be implemented by a computer program. Proceeding this way, we obtain analyt- 
ical expressions where D and n are arbitrary. Substituting numerical values into 
these formulas provides very precise numerical comparisons with other numerical 
methods. The results of the numerical methods at low n are very reliable and can 
be used to test the analytical results. If the analytical results pass this test, they 
can provide very accurate results at large n and, in turn, provide a test for the 
numerical stability of the numerical method. 

The simplicity of the calculations presented here relies on the fact that at the 
end, the arguments of the functions are set to zero. It will not be necessarily the 
case for higher order calculations and it is conceivable that some complicated sums 
will need to be performed numerically. Nevertheless, it is clear that the methods 
proposed here provide a drastic simplification when compared to a straightforward 
evaluation of the sums as they appear in Eqs. (3.6-9). 

In the infinite volume limit, one obtains more compact formulas. The coeffi- 
cients can be written as ratios of polynomials in c. The location of the poles at each 
order seems to obey some regularity. For the contributions which are calculable for 
arbitrary order, we see that when the order increases, some of the poles get closer 
to c = 2 and c = 4. This fact has to be kept in mind if one attempts to calculate 
the corrections to the gaussian approximation using a 1/D expansion. 
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APPENDIX 1: Calculation of the 6 n , m 

This is a Mathematica program which calculates the first four coefficients. 
The explanations between (* ...*) play no role . If this program is stored in a file 
named bees.math for instance, then it can introduced in a Mathematica session 
with the command << bees.math. When all the instructions are completed (it 
takes a few minutes), you will see completed. You can then use 63 for 63 etc... 
In order to check numerical values with the program given in Appendix 2, you 
can type for instance N[{bl, 62, 63, 64}/. {n — > 5, di — > 3}, 10], you should then 
obtain {1.075127131,0.8953433348,0.6726561277,0.4555210038}. This program is 
appended to the tex file. 

APPENDIX 2: Checking the Previous 
Calculation at Finite Volume 

This is a Mathematica program which calculates numerically the coefficients up 
to a given order for 2 n sites. It is convenient to take n < 8. The numerical values 
of n and D = di can be changed. The choice of values below, allows to check the 
numbers mentioned in Appendix 1. After entering the program in a Mathematica 
session, you get (after a few minutes) 1. + 1.0751271315 + 0.8953433348B 2 + 
0.67265612775 3 + 0.45552100385 4 . This program is appended to the tex file. 
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